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We are presenting here the new formulae for Bose-Einstein correlations (BEC) 

which contain effects of final state interactions (FSI) of both strong (in s-wave) 

and electromagnetic origin. We demonstrate the importance of FSI in BEC by 

^ , analysing data for e+e~ annihilation and for heavy collisions. The inclusion of 

vO 1 FSI results in the practical elimination (at least in e+e" data) of the so called 

^^ ' degree of coherence parameter A (which becomes equal unity) and the long range 

ff^ ' parameter 7 (which is now equal zero). 

(N 

O 

t^^ ' 1 Introduction 

r^ ' Bose-Einstein correlations (BEC) are the one of the most important current 

Mh, topics in high energy coUisions (in particular in heavy-ion reactions). One of 

(~^ ; the interesting problems present in the BEC is the physical meaning of the 

1^ , degree of coherence parameter A and the long range correlation parameter 7 

i-G ' which are usually introduced hy hand when analysing BEC data by means of 

the so-called standard formulald: 



N(±±)/N'''' = c[\ + \E2b]{1+iQ)- (1) 

The function E2B appearing here is called the exchange function and c is the 
so-called normalization parameter (also introduced by hand but we shall not 
discuss it here). 

Our approach to this problem is based on the observation that BEC for- 
mula Eq.(l), which has been obtained by using the plane wave approximation 
(asymptotic states) for both observed particles, should be corrected for the ef- 
fect of final state interactions (FSI) which can be of strong and electromagnetic 
type (Coulomb interactions). We have recently obtained several theoretical for- 
mulae for the BEC including effects of FSI of the strong type (in the s-wave, 
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isospin 1 = 2 channel) a, of the Coulomb intcractionsU'Q and of both types of 
FSI acting together B. 

The second order BEC without FSI is usually presented as the following 
convolution 

/■OO 

Jq 

/•OO 

= / p{R)p{r)\A,2\'d'Rd''r (2) 

Jo 

of the single particle source density functions p and the squared two particle 
amplitude 

Ai2 = —r=[exp{{ipi{xA - xi) + ip2{xB - X2)} 

+ exp{{ipi{xA -X2) +ip2{xB -xi)}] (3) 

which is symmetrized accordingly, cf. Fig. 1. 



Fig. 1 Identical boson exchange diagram. 
Assuming now the following Gaussian distribution for the source function, 
p{r) = e^^ /'^^ , one obtains Eq. (1) with E2B = e^^ "^ ^^ (and with c = 1, 
A = 1 and 7:^0). 

In the next section we shall present how to correct this formula for the 
presence of strong interactions FSI which are seen in the phase shift analysis 
of the TTTT and K^K^ correlations. The inclusion of Coulomb interactions, 
which are important when produced bosons are charged, demands however a 
profound change of the two particle amplitude Ai2. This is discussed in Section 
3. It is shown there also that, we can obtain information on the interaction 
range from data for unlike charged pairs tt+tt^. Section 4 deals with the most 
general case when strong and Coulomb FSI coexist together. Our concluding 
remarks are presented in the last section. 



2 Final state interaction in BEC for neutral particles 

FSI of the strong type are limited (due to the short range of strong interac- 
tions) to a small number of partial waves (in practice to s-wave only) and it is 
therefore sufficient to use the following expression for th^^mplitude describing 
the system of two identical bososns in their rest frame H'Q: 



J_p»(Pl+P2)-R- 
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V / „ikr I „ikr I „-ikr 



ikr 



(4) 



Here R = (ri + r2)/2, k — (pi — P2)/2, and r = (ri — r2) and 5 denotes t 
phase shift ^describing the corresponding FSI. The data of phase shifts of tttt 
and K^K^u, which were used in our analysis, are shown in Fig. 2. 



Fig. 2. Phase shifts for tttt and K^K^ collisions, (cf. [8], [9] for details.) 
After some algebra we obtain a following new formula of the BEC containing 
the strong FSI given by the respective phase shift function S:, 

^(2-)/^SG ^^^^ (e-^''3'/2Re[erfc(z)] + ^^^e~^''^'/'Re[evic{z)] 



8 sin S cos S 



/32Q2 
e-'^'Q'/^lm[eTic{z)]Y (5) 



where z = —if]Q/y/2, and where Re[erfc(— i/?(5/\/2)] = 1 relation was used 
(the degree of coherence parameter A was added by hand here). Notice that 
for (5 — > 0, i.e., when the FSI is switched off, we are recovering standard formula 
given by Eq. (1) with E = e~^ ^ /^ (if we add by hand additional parameters 
c and 7 in the form of long range correlation factor (1 + 7Q)). In Fig. 3 we 



present our analysis of BEC data for K^K^ pair production using the above 
approacha. 



Fig. 3. Analyses of the data of K^K^ pair. 

3 Final state interaction in BEC for charged particles - Coulomb 
interactions 



In the case of Coulomb type of FSI the amplitude A12 has to be described by 
the so called Coulomb wave functions: 
1 



A12 ^ 



V2 



[vl/(k,r) + *s(k,r)] 



(6) 



'i'(k,r) = r(l + i7y)e-^''/2e*''$(-i77;l;ifcr(l-cos6l)), 
*s(k,r) = r(l + i77)e-'^''/^e-*-''$(-i77;l;ifcr(l + cos6l)) 



where r — xi — 0:2, the parameter 77 — m,-cdC(/k, the momentum transLfer 
Q = {Pi ~ P2) — 2fc and $ denotes the confluent hypergeometric function $E£l. 
The BEC formula now reads: 
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(7) 
(8) 



where G{2k) = 27rr//(e^'^'' — 1) denotes Gamow factor, the first and the second 
parentheses in Eq. (g) correspond to the first and the second terms in Eq. (^ 



and (cf. ref. [3]) 

In{n,m) = An J dr r'+-+- p{r) , A„ = ^rUr\»V 

This formula, as it was shown injXcf.[3], can be also used to analyse data on 
unlike sign tt+tt^ pair productionEll. This data cover the Coulomb interaction 
region at small momentum transfer Q < AO MeV/c in P + Ta ^ tt^ +7r~ + X 
at proton cnerg}j_ZO GeV and have been so far analysed by using only Gamow 
factor correctioiO. However, when analysed by using instead our Coulomb 
formula with rj = mTra/2k, 

7v(+-)/7vSG^G(-ry)(l + Aic(-?7)) (9) 

they can provide us information on the interaction region in p + Ta — > 7r+ +7r^ 
+ X collisions which turns out to be about 2 fm, cf. Fig. 4. It should be noticed 
that the interaction range cannot be estimated from the analysis performed 
only by using the Gamow factor. 



Fig. 4. Analysis of data for the tt^tt" pair production. 
When one applies our formulae to single pion production in_.the external 
Coulomb field one founds an apparent quasi-scaling behavior E3 in the pion 
production yield N'^{(3, Z^ff) = 1 + Aic, namely one observes that 

N^iQ; /3, Ze//) « N^iQ; A x /3, A x Z^//), (1 < A < 3). (10) 

Because of this property, we have found that it is difficult to estimate the 
magnitude of nuclear fragments {Zf,ff) produced in the central region in heavy- 
ion collisions from the yield ratio A^'^ /N"^ only. 



4 Final state interaction including both strong and Coulomb in- 
teractions 

To describe a pair of the identical bosons including both strong and Coulomb 
interactions, we have to symmetrize the total wave function in the following 
way: 



1 



A,2 = ^[^ciKr) + n (k, r) + <^st (k, r) + ^^, (k, r)] . 



(11) 



Here ^c(k, r) denotes Coulomb wave function described above, superscript S 
denotes the symmetrization of the wave function and function <i>st (k, r) stands 
for the wave function induced by strong interactions. Assuming a source func- 
tion p{r) we obtain the following expression for the BEC in this case: 
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(12) 



Cst 



I.t = 
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{2ky^ exp (~i(77o + S^^^)) sin Jf V /r,2(1 + n)A2{n, 0) 



n=0 



/Ri(0)sin2 4'). 



Explicit expressions for quantities present here are given in ref.[5]. It is also 
shown there that our method is equivalent to the numericaJ solution of the 
Schodinger equation with strong and Coulomb potentials E3. Our final for- 
mula (containing three additional parameters: c, A and 7, added by hand) is 
therefore given as: 



Ar(±±)/ivBG(Q ^ 2k) 

X 1 + A, 



cfl 



Aic + Aec 
E2B 



1 



^IC 



iEC 



-'^Cs 



-^Cst + 1st) 

(1 + 7Q)- 



(13) 



It should be noted that the normalization c and an effective degree of coherence, 
i.e., the denominator of the ratio -E2b/(1 + Aic + AEc + -^Cst +-^st), are related 
to each other. For the sake of reference we shall also use in our analyses the 
standard formula as given be Eq. (1) (with E■i^^ exp(— /3^(5^/2)). We apply 
now our formulae to data for e+e^ annihilatioiESEj. Results of our analyses are 
shown in Fig. 5 and Table I. As seen in Table I, our (Eq.(13)) estimated values 
of the degree of coherence parameter A are systematically larger (approaching 



unity) than those obtained by the standard formula (Eq.(l)). Similarly, the 
long range correlation parameter 7 approaches now (approximately) zero. 



M 



Fig. 5. Example of the analysis (jLdata in e"^e annihilation for TPCt3 and 
OPALEJ collaborations. 
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ref.[14] 


0.92 ±0.06* 


0.61 ±0.05 


- 


- 


— 


Eq.(13) 


0.74 ±0.05 


1.10 ±0.04 


-0.00 ±0.02 


1.00 ±0.02 


44.2/35 


Eq.(l) 


0.91 ±0.06 


0.61 ±0.05 


0.08 ±0.03 


0.88 ±0.02 


41.0/35 


OPAL 












ref.[15] 


1.12 ±0.02* 


0.85 ±0.03 


- 


- 


336/73 


Eq.(13) 


1.09 ±0.04 


1.04 ±0.03 


0.00 ±0.00 


0.99 ±0.01 


124.4/74 


Eq.(l) 


1.34 ±0.04 


0.71 ±0.04 


0.04 ±0.00 


0.94 ±0.00 


118.7/74 



Table 1: Estimated parameters 



5 Concluding remarks 

We have presented several analytic formulae for BEG including the Coulombic 

and strong FSI obtained by us recently. Combining the seamless fitting method 
and the CERN MINUIT program in Eq.(13) we have analysed different sets 
of BEC data showing respectively: 

(1) The role of strong FSI by using data on K^K'^ phase shifts; 

(2) The ability to obtain by using our method the range of interaction from 
precise data on tt+tt^ production; 



(3) The possibility that A < 1 and 7 7^ values of parameters encountered in 
the standard analyses of data could be a reflection of the combined action 
of strong and Coulombic FSI which were not taken properly into account 
there. In fact, also the values of the source size parameter reported by 
various collaborations (after using the relation: f3 ~ v2R) and obtained 
by the standard formula Eq.(l) are systematically larger than values 
estimated by our formula Eq.(13). 
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Fig.2(a) 
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Fig. 2(b) 
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Fig. 4 
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